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Abstract 

We review the dual relationship between various compact U(l) lattice models and 
Abelian Higgs models, the latter being the disorder field theories of line-like topologi- 
cal excitations in the systems. We point out that the predicted first-order transitions 
in the Abelian Higgs models (Coleman- Weinberg mechanism) are, in three dimensions, 
in contradiction with direct numerical investigations in the compact U(l) formulation 
since these yield continuous transitions in the major part of the phase diagram. In 
four dimensions, there are indications from Monte Carlo data for a similar situation. 
Concentrating on the strong-coupling expansion in terms of geometrical objects, sur- 
faces or lines, with certain statistical weights, we present semi-quantitative arguments 
explaining the observed cross-over from first-order to continuous transitions by the 
balance between the lowest two weights ("2:1 ratio") of these geometrical objects. 
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1 Introduction 



In many physical systems, even the quahtative features along a line of phase tran- 
sitions may depend on the choice of the coupling constants which parametrize the 
bare interaction energies. Quite often one encounters the case where a continuous 
higher order transition changes into a discontinuous first-order one. In most cases, 
the cross-over behaviour is characterized by a tricritical point (TCP)|l|]. In field the- 
oretic language, the origin of a TCP is easy to understand on the basis of a Landau 
expansion of the effective potential in terms of a order- (or disorder-) parameter. 
One source for a TCP is the sign-change in the quartic term (assuming a stabilizing 
sextic term) in some range of the parameter space. A positive quartic term leads to 
a second-order transition, and a negative one to a first-order transition. A vanishing 
quartic term is associated with the TCP. Another important mechanism arises from 
the presence of a cubic term. Such a cubic term may be generated, in some range of 
the parameter space, by fluctuations of other fields. This always drives the transi- 
tion first-order. While the analysis of the Landau expansion is very elementary, the 
non-trivial problem is, of course, its derivation from the underlying physical system 
with quite complicated microscopic forces (the most prominent example being the 
Gorkov derivation ||^ of the effective Ginzburg-Landau theory of superconductiv- 
ity). It is therefore desirable to find alternative properties of the system which allow 
at least for a semi- quantitative understanding of the appearance of TCP's. In this 
lecture we exhibit such properties for compact U(l) lattice models, having in mind 
quite different physical systems (D=dimension): 

1. D=4 Quantumelectrodynamics (U(l)-LGT) 

2. D=3 Superfluid Helium (XY- model) 

3. D=3 Defect models of melting 

The common point of these systems is that they can be studied in 

a) the weak-coupling or "defect" expansion. This is formulated in terms of linesQ 
with long-range Biot-Savart-like interactions. This is why all three systems can 
be described, alternatively, by Abelian Higgs models in which the complex fields 
account for the line-like disorder of the systems and the gauge fields for the long- 
range interactions [0]. 

b) The strong-coupling or "stress" expansion. This is also formulated in terms of 
geometrical objects (surfaces, lines), but with no long-range interactions. This sim- 
plifies their statistical behaviour and will provide for the desired alternative criterion 
for the existence of TCP's. 

For the specific example of U(l) lattice gauge theory the different representations, 
being dual to each other, and their corresponding field-theoretic descriptions are 
summarized in fig.L The diagrams for the other systems are completely analogous. 



^When the dimension is reduced by one, D ^ D ~ 1, "lines" have to be replaced by "points". 



Figure 1: Interrelation between the dual "defect" (weak-coupling) and "stress" (strong- 
coupling) representations of U(l) lattice gauge theory and their equivalent field theoretic 
descriptions. 

2 The "Defect" Representation 

Let us first look at the weak-coupling expansion of the various U(l) lattice models. 
They all contain topological excitations with long-range interactions caused by the 
Nambu-Goldstone bosons of the U(l) symmetry. In a three-dimensional solid, the 
topological excitations are line-like defects, dislocations and disclinations IQ , and 
the Nambu-Goldstone bosons are elastic distortions (phonons), leading to a stress- 
field around each defect and to Biot-Savart-like interactions between defect elements. 
In superfluid Helium in three dimensions, the "defects" are vortex-lines and the 
long-range forces are caused by the superflow. In four-dimensional U(1)-LGT finally, 
the "defects" are world lines of magnetic monopoles P and the forces are due to 
electromagnetism. The upper right box in fig.l symbolizes the equivalent disorder 
field description of these "defect" lines in terms of a complex field \1/ interacting with 
a gauge field A^. In superfluid Helium and U(1)-LGT, the disorder field theory of 
line-like defects turns out 0, ^ to be an Abehan Higgs model. In a solid, it is a 
more complicated field theory of a similar type, involving two disorder fields, one 
for dislocation and one for disclination lines §]. 

Let us briefly recall the physical content of the dual equivalence. From Feynman's 
path-integral representation of quantum mechanics, we know that the statistical me- 
chanics of one fluctuating line ("orbit") corresponds to the quantum mechanics of 
one particle. It is then easy to see that the grand-canonical ensemble of fluctuating 
lines corresponds to the quantum mechanics of a many-particle system. This in turn 
is described most conveniently by a second quantized field theory. The long-range 
Biot-Savart-like nature of the interactions between line-elements is what permits 
their description by an Abelian gauge-field A^. The minimal coupling of to the 
disorder field \i/ leads then immediately to an Abelian Higgs model (scalar electro- 



dynamics). Indeed, in the lattice formulation of interacting defect lines, these steps 
can be carried out rigorously [|^ |[, leading to the Higgs action 



A 



-F. 



where = d^j^A^ — d^A^, and e,w?,g > can in principle be calculated from the 
couplings of the original compact U(l) lattice model, indicated by the circle in the 
center of fig.l. According to the famous Coleman- Weinberg argument |^, this 
Abelian Higgs model should always undergo a first-order transition when the mass 
parameter turns negative. The argument goes as follows: Assuming \1/ = const., the 
gauge fields can be integrated out, yielding a one-loop correction to ([l|) 



trlog(fc2 + e^|^n 
This amounts to additional interactions of the form 



1^1 



and 



in D 



(2) 



(3) 



1^1^ log 1^1 inD = 4 (4) 

In both cases, the additional term should drive the transition first-order. 

In three dimensions, this conclusion is in clear contradiction to the well estab- 
lished fact that the original D=3 XY model has a continuous phase transition 



111 , The crucial point in the whole argument is, of course, the assumption 



of an almost constant \l/-field. If this assumption breaks down, then also the gen- 
eration of a cubic term (in D=3) is no longer reliable and the transition may stay 
continuous even with A^-field fiuctuations. The decisive parameter is the ratio of 
the length-scales of A^- and fiuctuations. 



1 penetration depth(74^' 
a/2 coherence length 



(5) 



In three dimensions, it was possible to show that there exists a tricritical value of 
K ^ l/\/2, i.e. near the separation line between type-I and type-II superconductivity 
0. For small k, the Coleman- Weinberg mechanism is valid, leading to a first-order 
transition, whereas for large k, the transition stays second order. 

This led to the suggestion that also four- dimensional scalar electrodynamics could 



have a TCP [0. By the above duality arguments, this would also hold for U(l)- 
LGT. Indeed, the first-order nature predicted by the Coleman- Weinberg mechanism 
was at odds with early Monte Carlo investigations of the dual U(1)-LGT which all 
claimed evidence for a continuous transition and even reported estimates for critical 
indices [0. More recent work favors the existence of both transition regions, first- 
order and second-order, in agreement with ref.[13|. Only one paper claims the 



validity of the Coleman- Weinberg mechanism everywhere. The present status will 
be reviewed in more detail in the next section. 



The goal of this lecture is to present semi-quantitative explanations for these 
TCP's, working within the dual strong-coupling expansions of the various U(l) 
models. This has the advantage that, since the geometrical objects appearing in 
this expansion have no long-range interactions, subtleties of the type (|^)-(^) are 
absent. 



3 The "Stress" Representation 

Let us start by recalling the definitions of the various U(l) models mentioned in the 
introduction. The partition functions 

z= y: (n^) (6) 

{7— conf.} 

are products of local Boltzmann factors which may be chosen as 

r e^™'^® Wilson \ 

^ _ I g/3cose+7cos20 Mixed laction (7) 

i Ene-^^®-'""^' Villain J 

with standing short for 

r Vi7j - Vj7i = (Vi7j)^ Gauge 1 
Q = lVa XY \ model (8) 

[ Vi7i + V,7i = iyajf Melting J 

The product in (|^) runs over the plaquettes or links of the lattice, and Vj are the 
usual lattice derivatives (Vi/(x) = f{x + ?) — f{x)). The 7-configurations are, in 
U(1)-LGT, the euclidean electromagnetic fields, in the case of superfiuid Helium, the 
phase angles of the condensate, and in the melting model, the atomic displacements. 
The common important property of these models is the U(l) invariance 

7 ^ 7 + 27r (9) 

suggesting a Fourier ( "character" ) expansion of the Boltzmann factors 

00 

B{Q)= J2 Wte'"'' (10) 

6=— 00 

with Fourier coefficients 

Wb = r ^5(e)e-^^« (11) 

Inserting (|l^) into the partition function and summing over the 7-configurations, 
one obtains the strong-coupling expansion {Wi, = Wb/Wo) 

Z = ([[Wo) E (E^b) (12) 

{b— conf .} 



The algebraic structure of G in eq.(^ leads to conservation laws which constraint 
the admissible 6-configurations. It is most transparent to characterize them in terms 
of geometrical objects: 

{surfaces Gauge 
lines XY (13) 

complicated lines Melting 

In the U(1)-LGT, the 6-variables correspond to the electromagnetic field strengths, 
in the superfluid to the superfluid currents, and in the melting model to the physical 
stress. This is why we call the expansion (p!2D generically the "stress" -representation. 
The geometrical objects in ( [T3| ) are subject only to short-range contact interac- 
tions. Notice that the geometrical characterization of these "stress" -graph config- 
urations does not depend on the dimension D. This is in contrast to the "defect" - 
representation where the dimensionality of the geometrical objects depends on the 
dimension in which the duality transformation is performed. 



4 Analytical and Numerical Results 

Let us now analyze the "stress" representation (|12D in some detail. To each choice of 
action in (0) corresponds a special "natural" set of weights Wb in ([T2|). Their relative 
importance with increasing strength b = ±1, ±2, . . . can be studied conveniently by 
simulating each model with different actions and comparing their thermodynamic 
quantities such as their internal energies. 

We start with the comparison between Wilson's and Villain's action for which 
the weights are = Ib{P)/Io{P) {h '■ modified Bessel function) and = 

exp(— 6^/(2/?y)), respectively. The two actions are made as similar as possible by 
equating the b = ±1 weights which amounts to relating the Villain parameter /Jy to 
the Wilson parameter f3 as follows (see fig. 2) 

This relation was first written down by Villain ||15[ when analyzing the Wilson type 
of action of the XY model in terms of the discrete Gaussian model (Villain approx- 
imation). If furthermore the overall normalizations of the partition functions are 
adjusted, then the difference between the two actions lies all in the higher weights, 
W^iP) = h{(3)/W) + [/i(/3)//o(/3)]'' = exp(-6V(2/3y)) = for 6 > 2. In or- 
der to measure their importance, we have performed Monte Carlo simulations in the 
representation (^ with both Wilson's and Villains's action |T^. As a typical exam- 
ple, we compare in fig.3 our results for the internal energy of the U(l) lattice gauge 
model. The excellent agreement for low /3 up to the phase transition around /5 ~ 1 
demonstrates that, in this range, the systems are dominated by 6 = ±1 excitations. 
It is therefore not surprising that also the transition temperatures of the Villain 
action (with (5y as a free parameter) are mapped by ([T^ very precisely onto the 
corresponding ones of the Wilson action (see table 1). Furthermore, from fig.3 we 



Figure 2: The parameter (3v of the Villain approximation versus P of the corresponding 
Wilson action. The curve results from the requirement of equal weights for graphs of 
strength-1 in the "stress" representation for both actions. 



Figure 3: The internal energy of D=4 U(l) lattice gauge model with Wilson action in 
comparison with the Villain approximation. 



Model 


/3yc 


/5,(from VA) 


13c 


4D-U(1) 


0.643'* 


1.04 


l.OllP 


3D-XY 


0.33'= 


0.45 


0.4539^^ 


2D-XY 


0.73^= 


1.18 


1.12-^, 1.18f,1.14'* 



Table 1: Transition temperatures of U(l) lattice gauge model for D=4, and of XY 
model for D=2,3 compared with the values obtained from the Villain mapping (dy = 
— [21og(Ii(/3)//o(/5))]~"'^. The values are taken from: a) Jersak et al. [^], b) Gupta et 



al. |2§1, c) Dasgupta and Halperin d) Ferer et al. 0, e) Shugard et al. 0, f) 
Tobochnik and Chester g) Samuel and Yee [^], h) Harten and Suranyi 



W2/WI 


XY(/5 = 0.45) 


Gauge(/5 = 1.0) 


Wilson 
Villain 


0.11 
0.01 


0.24 
0.10 



Table 2: The "2:1 ratio" near the transition points of the D=3 XY and D=4 U(l) lattice 
gauge models. It shows that the "stress" graphs of strength-2 are much more important 
in the Wilson than in the Villain action. 



read off that at and above the phase transition the graphs of higher strengths prolif- 
erate much stronger using the Wilson action. Among these we expect the graphs of 
strength-2 to be most important. This suggests that the weight ratio W2/W1 should 
be the relevant distinguishing feature of the different actions. In the sequel, it will 
be called the "2:1 ratio". In the Wilson action, the "2:1 ratio" is much larger than 
in the Villain action as is shown in table 2. 

It is easy to convince ourselves that the higher ratios W^/Wi, . . . are practically 
irrelevant near the phase transition. For this purpose we compare the Villain action 
with further Monte Carlo simulations of the mixed /3 — 7 action (see eq.(|^)), which 
is chosen in such a way that the "2:1 ratio" coincides with that of the Villain action. 
The coincidence is achieved by requiring 

- expf-^) (15) 

, 4 

where Ib{P,'y) are generalized modified Bessel functions 

/,(/3,7) = r ^cos(60)e''^°'^(*)+^™^(2</,) ^^7) 

Here, it is convenient to solve eqs.(|l5|,|l6]) in the form P{Pv),l{l3v) which we shall 
call the "Villain locus" in the /3 — 7 plane. The comparison of our Monte Carlo data 



Figure 4: The internal energy of U(1)-LGT with Villain action. The lower data points 
are obtained from simulations of the f3 cos model transformed according to eq. (|T^ ) . The 
stars which fall practically on top of the Villain data are from a simulation with mixed 
action /3cos0 + 7 cos 20 treated according to eqs.( |l5tp^ ). 



for the internal energy of U(1)-LGT in fig.4 clearly demonstrates that the differences 
due to graphs of higher strength b = ±3, . . . (which still have different weights) are 
completely negligible - the data for the Villain action and the mixed /5 — 7 action 
(along the "Villain locus") fall almost perfectly on top of each other. 

As a conclusion, the differences between the Wilson and the Villain action near 
the phase transition are completely explained by the different weights of strength-2 
"stress" graphs, whose proliferation is much more pronounced in the Wilson case. 

That an admixture of strength-2 graphs with large enough weight can in principle 
drive the transition first-order, can be demonstrated easily in the disorder field 
theory of "stress" lines (i.e. in the mean-field formulation of the D=3 XY model 



with mixed action). This is described by an effective action ^1 



A = a|0ip + b\(f)i\'^ + c|02p - ^[4>l<Pt + c-c] + gradients + . . . (18) 

where the complex fields 0i, 02 represent strength-1 and strength-2 "stress" lines, re- 
spectively. The coupling 0i0i02 corresponds to the merging of two lines of strength- 
1 into one line of strength-2. The complex conjugate coupling describes the reversed 
process. Integrating out the field 02, one obtains an additional quartic term in 
01, — ^|0i|^, such that the total quartic term in 0i may change sign, signalizing a 
first-order transition. 

This explains at least qualitatively why in careful Monte Carlo studies with 
the Wilson action a first-order transition was reported while with the Villain action 



Figure 5: The parameters /?, 7 of the mixed action (3 cos + 7 cos 20 which can be studied 
by means of the improved Villain approximation ("Villain locus"). The fat and dashed 
lines show an interpolation to the critical points found by Jersak et al. j22U. The dotted 
line was estimated by those authors to be the locus of the Villain model in a different way 
from ours (which we believe to be less accurate). 



evidence for a continuous transition was claimed. Moreover, this picture is consistent 
with Monte Carlo studies of the mixed /? — 7 action by Jersak et al.[^ who located 
a TCP for slightly negative 7. Looking at fig.5, we see that our "Villain locus" 
indeed crosses the phase transition line in the range of second-order transitions. It 
would be interesting to investigate whether the TCP is connected with a universal 
tricritical "2:1 ratio" W2/W1. 

The most recent Monte Carlo renormalization group (MCRG) studies are in- 
consistent with this nice picture. They are interpreted as evidence for a first-order 
transition, with both the Wilson and the Villain action. A historic summary is 
compiled in table 3. The latest study [^] even speculates that the transition is 
always first-order, downto 7 = —00. The last result would imply that in the whole 
parameter range of the D=4 Abelian Higgs model which is covered by the dual /3 — 7 
action, the Coleman- Weinberg mechanism is indeed working. We feel, however, that 
even in view of this impressive list of work in table 3, the answer has not yet settled 
and much more high precision studies are necessary to provide the final answer. 

While in D=4 the situation is still quite controversial, in D=3 dimensions the nu- 
merical evidence is clearly against the analog of the Coleman- Weinberg mechanism. 



as advanced by Halperin,Lubensky and Ma ||10|- Here, it was confirmed many times 



Author 


7tcp 


Wilson 


Villain 


Bhanot (1982) gl[ 


> 


2 


2 


Jersak et al. (1985) ^ 


-0.11(5) 


1 


2 


Gr6sch et al. (1985) [gg 






metastability 


Morikawa et al. (1985) [|4] 


> 


2 


2 


Gupta et al. (1986) ^| 


> 


2 


2 


Burkitt (1986) |g 


> 


2 


2 


Lang (1986/87) |27| 


> 


2 




Lang and Rebbi (1987) |2§ 






2 


Decker et al. (1988) [|9l 




1 




A. Hasenfratz (1988) 


no 


1 


1 



Table 3: Summary of work on the order of the phase transition in the D=4 U(l) lattice 
gauge model with Wilson and Villain action and on the location of the tricritical point in 
the parameter space of the mixed /3 — 7 action. 



by different methods that both Wilson's and Villain's action undergo continuous 
transitions. The question was then whether there exists at all a parameter range 
in the mixed /5 — 7 action which shows first-order transitions. According to our 
above analysis, a good candidate was only the range 7 > 0. We therefore concen- 
trated on this range and found ||3l|, first in a mean-field (MF) treatment, first-order 
transitions for 0.166 < 7 < 0.375. The full MF phase diagram is displayed in 
fig. 6 (use the labelings on the right and top axis). Since already in the pure XY 
model (/3 = or 7 = axis), the MF transition temperatures are off by a factor 
0.45/0.33 = 1.36, we have rescaled the MF curves by this factor when comparing 
with our Monte Carlo simulations which we run in the range where the rescaled 
MF results show the largest entropy jump, 7 ^ 0.35. ..0.40. The numerical results 
confirmed the first-order nature of the transition. As a typical signal, we show in 
fig. 7 the double peak structure in the internal energy histogram, corresponding to 
tunnelings between two metastable states. Since the observed first-order transition 
is very weak and the 7-range is probably very small, we were not able to locate the 
two tricritical points with resonable accuracy. We can therefore only claim evidence 
for a short line of first-order transitions in the range 7 ~ 0.35. ..0.40. 

In defect models of melting, no TCP seems to exist. The reason is the very 
large activation energy of the graphs in the stress expansion so that there are no 
pretransitional excitations up to the point where the free energy intercepts the defect 
expansion. This is responsible for a jump in the slope of the free energy, i.e. for 
first-order transitions 1131. 



Figure 6: Phase diagram of the D=3 XY model with mixed action as obtained by mean- 
field methods and Monte Carlo simulations. 



Figure 7: Energy histogram near the first-order transition line of the D=3 XY model with 
mixed action. 
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